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The  Set  Packing  problem  has  a  dual  covering  problem,  which  asks  how  many  of  the  same  objects  are  required  to 
completely  cover  every  region  of  the  container,  where  the  objects  are  allowed  to  overlap.  Many  applications  arise  having 
the  packing  and  covering  structure.  Delivery  and  routing  problems,  scheduling  problems  and  location  problems,  switching 
theory,  wireless  network  design,  VLSI  circuits  and  line  balancing  often  take  on  a  set  covering  structure.  However,  if  one 
wishes  to  satisfy  as  much  demand  as  possible  without  creating  conflict,  it  takes  on  a  set  packing  format.  In  this  paper  a 
linearization  technique  is  developed  to  solve  set  packing  problems  with  linear  fractional  objective  function.  The 
correctness  of  the  algorithm  is  shown  by  an  numerical  example. 

AMS  Subject  Classifications:  90C10,  90C20. 

KEYWORDS:  Set  Packing  Problem,  Linear  Fractional  Set  Packing  Problem 
1.  INTRODUCTION 

The  Set  Packing  problems  are  the  interesting  topic  for  many  of  the  researchers  but,  surprisingly,  few  applications 
of  the  Set  Packing  formulation  have  been  reported  in  the  literature.  Some  of  them  are  described  below.  Ronnqvist  [19] 
worked  on  a  cutting  stock  problem  formulated  as  Set  Packing  Problem  and  solved  using  a  Lagrangian  relaxation  combined 
with  subgradeint  optimization.  Zwaneveld  et  al.  [24]  formulated  a  real  railway  feasibility  problem  as  Set  Packing  and 
solved  it  exactly  using  reduction  tests  and  a  Branch  &  Cut  method.  Kim  [12]  represented  a  ship  scheduling  problem  as  Set 
Packing  and  used  LINDO  software  to  solve  it.  Mingozzi  et  al.  [14]  used  an  SPP  formulation  to  calculate  the  bounds  for  a 
Resource  Constrained  Project  scheduling  Problem  using  a  greedy  method.  Rossi  [20]  considered  an  SPP  formulation  for  a 
ground  holding  problem  and  solved  it  exactly  with  a  Branch  &  Cut  method. 

The  railway  industry  is  rich  in  problem  that  can  be  modelled  and  solved  using  Set  Packing  format.  In  this  day  and 
age,  arguably  the  most  important  of  these  are  the  ones  that  concern  the  effective  allocation  and  utilization  of  available 
resources.  The  problem  of  routing  trains  through  railway  junctions  arises  at  each  of  the  level.  Railway  management  often 
face  the  task  of  deciding  between  a  number  of  possible  investment  alternatives  concerning  proposed  infrastructure 
modifications  to  junctions  of  which  the  most  influential  factor  in  making  the  final  decision  is  capacity.  Railway 
management  are  very  interested  in  knowing,  with  precision,  what  level  of  rail  traffic  the  modified  infrastructure  would 
cater  for.  This  effectively  involves  determining  the  maximum  number  of  trains  that  could  be  routed  through  the  junction 
within  a  given  time  horizon.  Scheduling  airline  flight  crews  to  airplanes  is  another  application  of  set  packing.  Each  airplane 
in  the  fleet  needs  to  have  a  crew  assigned  to  it,  consisting  of  a  pilot,  copilot,  and  navigator.  There  are  constraints  on  the 
composition  of  possible  crews,  based  on  their  training  to  fly  different  types  of  aircraft,  as  well  as  any  personality  conflicts. 
Given  all  possible  crew  and  plane  combinations,  each  represented  by  a  subset  of  items,  we  need  an  assignment  such  that 
each  plane  and  each  person  is  in  exactly  one  chosen  combination.  After  all,  the  same  person  cannot  be  on  two  different 
planes,  and  every  plane  needs  a  crew.  We  need  a  perfect  packing  given  the  subset  constraints. 
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There  are  many  applications  of  fractional  packing  problems  in  real  life  as  well.  Some  of  them  are  'air  line  crew 
scheduling',  'truck  routing',  'political  districting',  'information  retrieval',  etc.  For  example,  suppose  an  air  line  company  has 

'  m  '  flights,  /  =  jl,  2, ...  ,  m}  to  operate  upon  and  '  n  '  crews,  J  =  jl,  2, ...  ,  ft}  at  its  disposal,  it  being  understood  that 
a  crew  can  handle  at  least  one  flight.  Let  Cj  >  0  be  the  profit  earned  by  the  company  when  its  j'h  crew  is  in  operation 

and  let  d  .  be  the  some  utility  of  the  jth  crew  paid  by  the  company.  Also  a  >  0  be  the  fixed  value  of  this  utility  function. 

Now  the  company  is  interested  in  scheduling  its  crew  in  such  a  way  that  all  the  flights  are  covered  and  the  ratio  of  the 
earning  and  the  utility  function  is  maximized. 

If  a  variable  X  ■  is  defined  as 

fl  if  fh  crew  is  in  schedule  Jl    if  ith  flight  by  j"1  crew 

Xj  —  s  and  djj  =<  ,  then  fractional  set  packing 

[0  otherwise  [0  otherwise 

problem  is  to  find  a  set  of  crews  /  <Z  /  that  covers  all  the  flights  and  Maximize   ;  je/, 

YsdjXj+a 

2.  THEORETICAL  DEVELOPMENT 
2.1  Set  Packing  Problems 

Consider  a  set  /  =  {1, 2, m]  and  set  P  =  {Pl ,  P2, Pn  }  where  P.  C  /  and  j  e  /  =  {1, 2,...,  n\ .  A  subset 
J  of  J  is  said  to  be  a  pack  of  /  if  |J  Pj  =  I, and  j&k  e  J* ,  j*k  and/>  f]Pk  =  </>■  Let  a  weight  C  .  >  0  be 

associated  with  every  je/.  The  total  weight  of  the  packing  J  is  equal  to  ^  C  ■. 

The  linear  set  packing  problem  (PkP)  is  to  find  a  pack  of  maximum  weight  subject  to  the  condition  that  at  the 
most  one  of  the  utility  is  satisfied.  Mathematically  the  problem  is 


(PkP)  Max  f{x)  =^JCjxj 


n 

Subject  to  X  d-ijXj  <l,  iel  (2.1) 

7=1 

Xj  =  0  or  1 ,  j  gJ  (2.2) 

fl         if  j  is  in  the  pack  fl         if  i  e  Pj 

Where    x ,  =  <  and    aif  =  < 

[0         otherwise  [0  otherwise 

Analytically,  the  Linear  Fractional  Set  Packing  Problem  (LFP)  with  the  same  restrictions  is  to  find  a  pack  of 
maximum  weight:  therefore,  mathematically,  the  problem  is- 
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(LFP)  Max  /(*)=  — 

j=i 

n 

Subject  to  X  a.,jc,  <  1,  i  e  /  (2.3) 

7=1 

Xj  =0  or  1,  j  gJ  (2.4) 

fl         if  j  is  in  the  packing  f  1         (/"  i  e  P 

Where    x .  =  <  and  a..  =  < 

[0         otherwise  [0  otherwise 

n 

It  is  assumed  that  C}  '  S  and      '5  are  non-negative  and  a  is  a  scalar  such  that  ^idjXj  +  a  >  0  .  In  matrix 
form,  (LFP)  can  be  written  as 

Maximize  f(x)  = 


Dx  +  a  ■ 
subject  to  Ax  <b 

Where      xT  —  (jc, ,  JC2,...,  JCn)  with  X-  =  0  or  1,  j  =  1, 2,..., n.     Here      C  =  (c,,c2,...,c„)  e  7?"  and 

D  =  (dl,d2,...,dn)  e  i?"  are  row  vectors,  A  is  an  mxn  matrix  of  zeros  and  ones  and  bT  =  (1,1,..., 1)  is  a  row  vector  of 
ones. 

DEFINITIONS 

•  Pack  Solution:  A  solution  X  which  satisfies  (2.1)  and  (2.2)  is  said  to  be  a  pack  solution. 

•  Redundant  Pack:  For  any  pack  J ,  a  column  j   G  J  is  said  to  be  redundant  if  /  +  j  j  J    is  also  a  pack. 

•  If  a  pack  contains  one  or  more  redundant  columns,  it  is  called  a  redundant  pack. 

•  Prime  Pack  Solution:  A  pack  J    is  said  to  be  a  prime  pack,  if  none  of  the  columns  corresponding  to  J  e  J  is 
redundant.  A  solution  corresponding  to  the  prime  pack  is  called  a  prime  pack  solution. 

•  Pseudo  concave  Function:  Let  /be  a  differentiable  function  defined  on  an  open  set  T  CZ  Rn .  Let  S  CZ  T 

and  Xl,X2  e  S ,  then  /  is  said  to  be  pseudo  concave  if  Vf(X2)T(Xl  —  X2)  <  0  =>  /(Xj)<  f(X2). 

2.1.2  Following  are  the  Results  that  Form  the  Basis  of  the  Algorithm  to  Enumerate  the  Problem 

Theorem  l.[5]  If  J   =  { j  :  X} ■  =  1}  is  any  prime  pack  of  (LFP)  then  X  =  {x} }  is  an  extreme  point  of  the  convex 
set  formed  by  feasible  region. 

Theorem  2.  If  the  objective  function  in  (LFP)  has  finite  value  then,  there  exists  a  prime  pack  solution  where  this 
value  is  attained. 
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Proof:  Let  a  finite  optimal  solution  of  (LFP)  exists  at  XQ  e  S  then  the  optimal  value  is 
Cx 

/(X())  =  _C^L_. 
DxQ  +  a 

Let  Jo  be  the  pack  corresponding  to  the  solution  X0  .If  Jo  is  the  prime  pack,  then  it  is  done,  otherwise  a  prime 

pack  can  be  derived  from  Jo  by  considering  the  redundant  columns.  Let  /[  be  the  prime  pack  obtained  from  Jo  and  Xl 
be  the  corresponding  solution  of  (LFP)  such  that 

Dxx  +  a 

Since  C, D  >  0,  a  is  the  positive  and  J{  ID  Jo,  therefore, 

> 


Dxx  +  a     DxQ  +  a 
Or  f(Xl)>f(x0) 

As  f(x0)  is  the  optimal  value  of/(x)  ,  therefore,  <  f(x0)  .  Hence  f(xY)  =  f(xQ)  .  Which  proves 

that  there  exists  a  prime  pack  solution  yielding  the  optimal  value  of  the  objective  function  of  (LFP). 

Theorem  3.  Let  f(x)  be  a  pseudo  concave  function  defined  on  feasible  set  S  and  X  eS  then  X  is  an  optimal 
solution  for  the  program 

Maximize  f(x) 

* 

if  and  only  if  X   is  an  optimal  solution  for  the  program 

Maximize  Vf(x*)Tx 

where  S  is  the  feasible  region. 

Proof:  Let  X  be  an  optimal  solution  for  the  program  (LFP),  therefore,  f(x)>  f(x) ,  V  xeS. 

As /is  differentiable  at  X  ,  therefore, 
f{x)  =  f{x*  +  x-x*) 

=  f(x*)  +  Vf(x*)T(x  —  x*)  +  a(x*  ,x-  x*)\x  —  jc*| 

Where    a(x* ,  X  -  x  )  — >  0  as  X  — >  X* 

Since     f(x*)  >  f(x)  =>  Vf(x  )T(x  -  x*)  +  a(x* ,  x  -  x*) \x  -  x*\ <  0 
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And       a{x* ,  X  -  X*)  — »  0 

^Vf(x*f(x-x*)<0 

^Vf(x)Tx<Vf(x)Tx 

=>  X  is  an  optimal  solution  for  the  program  MaxVf(x  )T  X  . 

Conversely,  let  X  be  an  optimal  solution  of  MaxVf(x*)T  X  therefore,  Vf(x*)TX*  >  Vf(x*)  X,VxeS 
^V/(x)(x-x)<0,Vxg5 

Since  f  is  a  pseudo  concave  function,  therefore,  f(x)<  f(x  ),VjceS. 
Hence  X  is  an  optimal  solution  for  (LFP). 

Following  is  the  algorithm  developed  to  enumerate  the  given  set  packing  problem. 
ALGORITHM 

Step  1:  Consider  a  Linear  Fractional  Set  Packing  Problem  (LFP).  Form  the  corresponding  continuous  program 

(LFP')  by  embedding  the  feasible  region  into  R"  (a  cube  with  n  vertices).  Let  S  be  the  feasible  set  for 
(LFP'). 

Step  2:  Choose  a  feasible  solution  XQ  e  S  such  that  V/(x0)  ^0.  Form  the  corresponding  linear  program 
(LP). On  solving  (LP)  let  Xjbe  its  optimal  solution.  If  X1  =  XQ  then  this  is  the  required  solution  of  the 
given  problem,  otherwise  let  S1  ={xl}. 

Step  3:  Starting  with  the  point  X1 ,  form  corresponding  (LP),  let  its  optimal  solution  be  X2  ^  Xl  .  Update  Sl  i.e. 
^l  =  l-^i'^}- 

Step  4:  Repeat  step  3  for  the  point  X2  and  suppose  at  the  ith  stage  St  =  {xl,X2,...,Xi }.  Stop  if  at  the  (i+l)th  stage 
XM  G  Si ,  then  XM  is  the  optimal  solution  of  (LFP). 

Step  5:  If  XM  is  an  optimal  solution  of  the  form  0-1  then  it  is  a  solution  of  (LFP)  otherwise,  go  to  Step-6. 

Step  6:  Apply  Gomory  cuts  to  find  a  solution  of  the  0-1  form  and  the  corresponding  prime  cover. 

Note:  The  algorithm  must  terminate  after  finite  number  of  steps  as  it  moves  only  on  the  vertices  of  the  feasible 

cube-/?" ,  which  are  finite  in  numbers,  i.e.  convergence  is  must. 
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NUMERICAL  EXAMPLE 


Max  f(x) 


(LFP) 


2x{  +  9x2  +  6x3 
x1  +  4x2  +  x3  +  2 
subject  to        xl  +  x2  <  1 
x2  +  x3  <  1 


xx  +  x3  <  1 
x1,x2,x3  =  0  or  1 

where      /=  {1,2,3}  ,  /  =  {1,2,3} 
Step  1:  The  corresponding  (LFP')  is 

rt        2xl  +  9x2  +  6x3 
Max  f(x)  = — 1  

x  =  (x, ,  x2 ,  x3 )  e  S  =  {(Xj ,  x2 ,  x3 )  I  x1  +  x2  <  1,  x2  +  x3  <  1,  xt  +  x3  <  1,  xt ,  x2 ,  x3  >  0  } 

Step  2:  Choose  X0  =  (1,0,0)  as  one  of  the  feasible  solution  of  (LFP')  with  y/"(x0)  ^  0 
The  corresponding  (LP)  is 

™    ■   ■  \T      4       19  16 

Maximize  V/  (x0 )  x  =  —  x,  +  —  x2  +  —  x3  :  x  e  6 

Now  it  is  a  linear  problem  therefore  can  be  solved  by  simplex  method  i.e. 
After  applying  the  simplex  algorithm  the  final  optimal  table  is  as  follows: 

Table  1 


4/9 

19/9 

16/9 

0 

0 

0 

CB 

B 

xB 

Yi 

Y2 

Y3 

Si 

s2 

S3 

19/9 

x2 

Vi 

0 

1 

0 

Vi 

Vi 

-1/2 

16/9 

x3 

Vi 

0 

0 

1 

-111 

Vi 

Vi 

4/9 

Xi 

Vi 

1 

0 

0 

Vi 

-111 

Vi 

39/18 

0 

0 

0 

7/18 

31/18 

1/18 

Hence  the  optimal  solution  is  (1  /  2, 1  /  2, 1  /  2)  ,  which  is  not  of  the  form  0  or  1,  therefore  apply  the  Gomory  cut 
to  get  integer  solution. 

-11         1  1 

The  cut  is  S4=  h  —  X,  +  —  Xs  +  —  X, 

2     2       2  2 

Now  inserting  this  additional  constraint  in  the  optimal  simplex  table  the  next  iterative  table  is 

Table  2 


c. 

4/9 

19/9 

16/9 

S3 

s4 

CB 

B 

xB 

Yi 

Y2 

Y3 

Si 

s2 

19/9 

x2 

Vi 

0 

1 

0 

1/2 

Vi 

-1/2 

0 

16/9 

x3 

1/2 

0 

0 

1 

-1/2 

Vi 

Vi 

0 

4/9 

Xi 

Vi 

1 

0 

0 

1/2 

-111 

Vi 

0 

0 

s4 

-1/2 

0 

0 

0 

-1/2 

-1/2 

-111 

1 

39/18 

0 

0 

0 

7/18 

31/18 

1/18 

0 

Now  apply  the  dual  simplex  method,  drop  S4  and  enter  S3  next  iterative  table  is 
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Table  3 


4/9 

19/9 

16/9 

0 

0 

0 

0 

CB 

B 

xB 

Yj 

Y2 

Y3 

Si 

s2 

s3 

s4 

19/9 

Y2 

1 

0 

1 

0 

1 

1 

0 

-1 

16/9 

Y3 

0 

0 

0 

1 

-1 

0 

0 

1 

4/9 

Yj 

0 

1 

0 

0 

0 

-1 

0 

1 

0 

s3 

1 

0 

0 

0 

1 

1 

1 

-2 

19/9 

0 

0 

0 

1/3 

5/3 

0 

1/9 

Hence  the  optimal  solution  of  (LP)  is  X1  =  (0,1,0)  ,  which  is  not  equal  to  X0  =  (1,0,0)  therefore  let  Si  =  {  Xl  }. 
Step  3:  Now  starting  with  the  point  X{  the  corresponding  (LP)  is 

Maximize  Vf(x,)Tx-—x,  +  —  x,  +  —  x,  :  x  e  S 
1         36  1    36  2    36  3 

Now  it  is  a  linear  problem  therefore  can  be  solved  by  simplex  method. 

After  applying  the  simplex  algorithm  the  final  optimal  table  is  as  follow: 


Table  4 


Cj 

3/36 

18/36 

27/36 

0 

0 

0 

CB 

B 

XB 

Yi 

Y2 

Y3 

Si 

s2 

s3 

0 

Si 

1 

0 

2 

0 

1 

1 

-1 

27/36 

x3 

1 

0 

1 

1 

0 

1 

0 

3/36 

0 

1 

-1 

0 

0 

-1 

1 

27/36 

0 

6/36 

0 

0 

24/36 

3/36 

Hence  the  optimal  solution  of  (LP)  is  X2  =  (0, 0, 1)  ,  which  is  not  equal  to  Xl  =  (0, 1, 0) ,  therefore,  let 


Now  starting  with  the  point  X2  the  corresponding  (LP)  is 


Maximize  Vf(x2)Tx-—x2  +  —  x3  :  x  e  S 


Now  it  is  a  linear  problem  therefore  can  be  solved  by  simplex  method. 


After  applying  the  simplex  algorithm  the  final  optimal  table  is  as  follows: 


Table  5 


C| 

0 

3/9 

12/9 

0 

0 

0 

CB 

B 

XB 

Yi 

Y2 

Y3 

Si 

s2 

s3 

0 

Si 

1 

1 

1 

0 

1 

0 

0 

12/9 

x3 

1 

0 

1 

1 

0 

1 

0 

0 

s3 

0 

1 

-1 

0 

0 

-1 

1 

12/9 

0 

1 

0 

0 

12/9 

0 

Hence  the  optimal  solution  of  (LP)  is  X3  =  (0,0,1) ,  which  is  equal  to  X2 
solution  for  the  original  (LFP)  with  optimal  value  2. 

CONCLUDING  REMARKS 


—  (0, 0, 1)  therefore  this  is  the  optimal 
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In  this  paper  we  have  considered  the  Set  Packing  Problem  with  linear  fractional  objective  function  the  technique 
developed  is  enumerative  one  which  go  through  the  vertices  of  the  given  problem  and  gives  the  optimal  solution.  Most 
importantly  the  algorithm  just  needs  the  pseudo -concavity  of  the  objective  function. 
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